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It is shown that the vacuum of QED2 in Minkowski spacetime does 

£""\ I not favour a periodic electric mean field. The projected effective action 

exhibiting a genuine dependence on the non-vanishing background field 

*vi has been introduced. The functional dependence of the energy den- 

sity of the vacuum on the assumed periodic vacuum expectation value 
of the vector potential is determined from the component T 00 of the 
energy-momentum tensor at one- loop order. Treating the background 

^^ \ field non-perturbatively, the energy of the vacuum in the presence of a 

C^ ' periodic mean field is found not be equal to the negative of the effective 

^—i ■ action. 
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It is well-known that the perturbative vacuum of QED is trivial, i.e. it is 
characterized by the vanishing expectation value of the electromagnetic vec- 



tor potential. Our goal is to investigate the vacuum non-perturbatively in the 
framework of the continuum theory in Minkowski spacetime. As a working hy- 
pothesis we assume the presence of a periodic electric background field. The 
interaction of the fermions with the background field is treated exactly by inte- 
grating out the fermion fields, whereas the quantum fluctuations of the electro- 
magnetic field are taken into account at one-loop order. The non-perturbative 
treatment of the interaction with the background introduces infinitely many 
non-renormalizable, i.e. irrelevant interaction vertices. It may happen that 
such vertices generate non-trivial, periodic (otherwise anti-ferromagnetic) vac- 
uum structure. Similar examples are known for Yang-Mills theories on the 
lattice 0, where the existence of various anti-ferromagnetic vacua have been 
established due to irrelevant interaction terms. 

Our problem setting can also be motivated by the following intuitive pic- 
ture. In the presence of a background scalar potential A = acos(£x) with 



overcritical amplitude a > m (with the rest mass of the electron m) local- 
ized electrons of the Dirac sea may tunnel into localized electron states of the 
same energy above the mass gap (see Fig. [I]). Therefore, the creation of a 
certain amount of electron-positron pairs with the spatial separation of n/£ is 
imaginable. This might stabilize itself with a periodic charge density (but zero 
net charge) and a periodic mean electric field. The expected effect might be 
destroyed due to the uncertainty principle. Localisation in an interval U/2£ 
means a momentum and, consequently, an energy spread of order AE = j| Ap 
and the effect may occur only if a >> m + 



AE 
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Figure 1: Intuitive picture of pair creation in overcritical periodic external 
electric field 



2 Projection method 

It is not completely trivial to define an effective action depending on a back- 
ground field A^{x). The simple shift of the integration variable A^(x) to 
a^{x) = A^{x) - A^{x) is not the answer requested, since a dependence on 
the background A^(x) shall only occur due to terminating the loop expansion. 



Below we introduce an effective action exhibiting a genuine dependence on the 
background field. 

Our method of defining the sector of QED for field configurations belonging 
to quantum fluctuations around a given vacuum expectation value (A^(x)) 
consists of the following steps: 

1. Projection to the sector with a given background field. The generating 
functional of QED in Minkowski spacetime has the form: 

Zqed = Jv^V^VAexp{iS em [A^}+iS D [A,^,^}} (1) 

with the vector potential A^(x) (fi = 0, 1, 2, 3) and the electron-positron 
field ip(x), where S^fA, £] and Sd[A, V','0] are the action of the electro- 
magnetic field and the Dirac action, resp., in the covariant gauge with 
the gauge parameter £. Let the vector n M (x) be introduced in the space 
of the vector potential configurations. Multiply the integrand of the 
generating functional by the factor 1 written in the form 

1 = f dc S( j dxA^n" - cQ) (2) 

with Q = TV the spacetime volume. Then, we find 

Zqed = J dcZ' QED [n,c], (3) 

where the functional 

Z'qedI^c] = Jv$VipVAexp{iS em [A,Z]+iS D [A,4>,ij;]} 

5( f dxA^n" - eft) (4) 

can be concieved as the generating functional of the sector belonging to 
vector potential configurations in a hypersurface orthogonal to n^(x). 
The projected generating functional Z l Q ED [n 1 c] is gauge invariant for the 
choice of n^(x) satisfying the condition d ll n^{x) = 0. The projected 
effective action S e ff [n, c] is defined by 

ZqedI^c] = exp{LS e// [n,c]}. (5) 

For the choice n^{x) = A^(x) with a given background field A^(x) satis- 
fying the Lorentz condition d ll A IJ '(x) = 0, we obtain the effective action 



for the sector belonging to the given background A^(x). Introducing 
the shifted integration variable a M (x) = A^(x) — A^{x), the projected 
generating functional can be rewritten as 

z 'qed\^A = J dX J Vaexp{iS em [A + a,£}} 

exp jiA ( i dxA^Af, + I dxA^a^ - cfij 1 
Jv^jVtljexp{iS D [A + a^^}}. (6) 

2. Identification of the background with the vacuum expectation value of the 
vector potential. As to the next, it is required that 

(A^(x)) = A^{x), i.e. (a") = 0. (7) 

This condition is used to determine the constant c as the functional of 
the vacuum expectation value of the vector potential (^(x)). 

For later use it is useful to introduce the external sources j M (x), fj(x), and 
r](x) coupled to the quantum fluctuation a^{x) of the vector potential, and to 
the fermion fields ip(x) and ip(x), resp. Then, we find instead of Eq. @ the 
expression 

Z QED[A,c,j,r],fj] = J dXj VaexpliS\[A + a,£,c] + J dxfa^ 

Zf[A,v,v} (8) 

with 

S x [A + a,\,£,c] = S em [A + a, £] + A ( f dxA^A^ + f dxA^a^ - cQj ,(9) 

and 

Z F [A,r),fj] = Jv^Vijexp{iS D [A + a,^,ij]}. (10) 

3 Projected effective action 

The projected generating functional (|6|) defines the projected effective action 
via Eq. (|5|). We determine it in one- loop approximation, i.e. we replace the full 
vector potential A^(x) in the Dirac action Sd by the background field A^(x). 



Thus, the generating functional is factored into an electromagnetic part and a 
fermionic part, 

Z'qed[A> c >3>V,v\ = z 'em[A,c,j}Z F [A,r),fj}. (11) 

The explicit forms of the actions are given as 

S ern [A + a,$ = ^{AD- l A)+ l -{aD- l a) + {aD- l A), (12) 

and 

S D [A] = (i'G- 1 ^) (13) 

with the inverse of the photon propagator in Lorentz gauge, 

D^(x,y) = [g, v U x +(^-l)dld:]5{x-y), (14) 

and the inverse of the fermion propagator in the background field A^(x), 

G-\x,y) = (i 7 ^-a M (x))-m)5(x-y) (15) 

with the Dirac matrices 7^ (/i = 0,1,2,3). For the sake of simplicity, the 
notation (fOg) = J dxdyf a (x)O a b(x,y)g b (y) is used, where a and b are either 
Lorentz or spinor summation indices. In the one-loop approximation the path 
integrals are Gaussian ones and can be performed explicitly, leading to 

1 1 

\nZ' em [A,c,j] = -- Tr lnZT 1 - -ln( AD A) + '-(jD.j) - (jA) 

+ 2jADA) +Cn jADAy (16) 

and 

\nZ F [A,rj,fj] = Tr InG' 1 + i(fjGr]) (17) 

with the modified photon propagator 

DT(x,v) = D~(x,y) - I^MU^ IdvD-i^Mv) (lg) 

Now we restrict our considerations to 1 + 1 dimensional systems and time 
independent periodic backgrounds of the form 

A tl = 5^acos{£x 1 ) (19) 



satisfying the Lorentz condition. 

The constant c can be determined by using the fact that 

= /Dot dX— exp | iS\[A + a, £, c] + / dxfa^ \ . 

From this we find that 

c Q = / dxA^. 

One establishes now for the one-loop effective action: 

iS eff = lnZ f em [A,c ,j = 0]+\nZ F [A,r ] = 0,r ] = 0] 

1 1 2 

= Tr InG" 1 -- Tr InD' 1 + n\a 2 £ 2 - -\n^-. 

2 4 2 t 2 



In the infinite volume limit one finds 



k k z ^ 



(20) 
(21) 



.'22) 



(23) 



6 

4 


2 5 H H S 
it j: is it i 
ii Is s t Is i- 
s t t t s t t t t 
tit t t Si \ i 
* 1 1 t i it 1 1 
: tt tt :: :: 

i i 1 1 1 

£222* 

It It It It It 
t Is 1 1 IS IS 

t t t f 1 1! i I 

ifslSiSfS 
SI t t 1 1 1 1 t. 
tt ts tt ts t 
ts ts tt tt t 
i t s s 1 






w o 








-2 


h •••••••""""**" "***"™ 4 ™ilImIS 


-4 
-6 


•••••«:::n«"n!«!*"*" "***!!!in""n;«at"«»» 

***^]**^";;s::tt»»«»»»»»tt:K:;;;n«!i""** 
••••wt::;s«^«!"* * ""n*"*«"** ,,,M ** 



-0.6 -0.4 



-0.2 0.0 

k 



0.2 0.4 



Figure 2: Energy eigenvalues with changing momentum 



Here the first and the second terms represent the energy density of the 
Dirac vacuum and that of the free electromagnetic field, resp. (see e.g. 0). 
J2 > denotes the summation over all non-negative single fermion energies e& > 0, 
being the energy eigenvalues of the Dirac equation in the external field A^{x). 
The third term on the r.h.s. of Eq. (^) is just the negative of the energy 
density of the periodic electric background field. The last term of the effective 
action fl22[) gives a vanishing contribution to the action density in the infinite 
volume limit. 

The meaning of the first two terms of Eq. (^3|) might lead one to the false 
conclusion that the negative of the action density is equal to the energy density 
of the vacuum, as it would happen if the background field were constant J| . 
Determining the energy density from the energy-momentum tensor we will 
show that the energy density does not equal the negative of the action density 
for inhomogeneous background field. 




Figure 3: The difference of the negative of the effective actions in the presence 
and in the absence of a periodic background field for various values of a and £. 



As 1 + 1 dimensional QED is a superrenormalizable theory, the action den- 
sity is UV finite after subtracting the action density of the free vacuum S e ffj ree , 
i.e. that of the vacuum in the absence of the background field. This difference 



— r = fi _1 (— S e ff + S e ffj ree ) has been calculated numerically by solving the 
Dirac equation for the single-fermion energies. 

These energy eigenvalues show the same band structure plotted against the 
momentum in fig. |2| as the electrons in the Kronig-Penney model || because 
of the periodic potential. Fig. |3| shows how — T changes by modifying the 
amplitude a and the wavelength £ of the potential. (The fermion rest mass is 
chosen for m—1.) One can recapitulate from fig. § that the surface —T(a,£) 
has only a single stationary point at the origin of the parameter space (a,£), 
i.e. the path integral defining the vacuum-vacuum transition amplitude is 
dominated by the trivial, identically vanishing field configuration A^(x) = 0. 



4 Mean field equation 



The vacuum expectation value of the current can be written at one-loop order 
in the following form: 



{^(x)^ip{x)) 



r/A / Vae iSx 
S 



S 



iSA^x) 



e ~' lS eff 



iSA^x) 



js eff _ z 



dX / Vae lbx Z F 



S 



F~ 



i5AK( x ) Zem 



(24) 



The first term on the r.h.s. of Eq. fl24|) is just the first functional derivative 
of the effective action 5S e ff/5A^{x). In the infinite volume limit the second 
term gives: 



1'Ti 



s 



exp 



[ 2 (ADA) J iSA^x 
A u (x) ( 5c 



exp 



2(ADA) J 



(DAY(x) ySA^x) 
£ 2 A»(x) = d v F^ 



- A"(x) 



(25) 



with the field strength tensor F^ u evaluated from the background A 11 . Here 
we used that 



5c 



SA^(x) 



2#(i) 



due to Eq. fl2"U|). Thus, one obtains the equation 

5S eff 



{^(x^i'ix)) 



SA^x) 



+ d v F 



)J,V 



(26) 



(27) 



For the background field configuration making the effective action extremum, 
one recovers the vacuum expectation value of the Maxwell equation. Thus, 
the effective action has an extremum at the background field configuration 
coinciding with the mean field solution. In our numerical search, this is the 
trivial extremum found at a = £ = 0. 

Satisfying the necessary condition of the extremum of the effective action, 
Eq. ( p7|) results in the Poisson's equation for the mean field A , which must 
be considered together with coupled set of operator equations for the quantum 
fields ip(x), ip(x) and a M (x). The latter equations should be solved and the 
result substituted in the l.h.s. of Eq. fl2~T|), in order to make the charge density 
explicit. 

5 Energy density of the vacuum 

The symmetric energy momentum tensor is defined by ||: 
dC dC dC 

with the Lagrange density 

c = -\f^f^ - hd^f + XA ^ - Ac o 

+ l - (^D^ - D$-fil>) - m^ (29) 

corresponding to the action / dxC = S\ + Sp and with the covariant derivative 
D» = d^ - iA». 

Substituting the ansatz ([H^) one obtains for the component T 00 of the 
energy-momentum tensor: 

T 00 = TZ,M+TZAAa,X}+TZ,o[A}+T x °[A,c] 

+Tf[A + a,$,iJ;] (30) 

where TZ a ( a = 0, 1, 2) are the terms independent of the fermion field and 
being of the order (a) a . Due to the constraint (|7|), the expectation value of 
the first order term vanishes, therefore it can be neglected. The second order 
term 

TZM = -d a fl d°a^ + 2d a fl d^a -d fl a d^a 

+^d fi a v d' J 'a u - ^a^af + ^-d^d v a v (31) 



is the expression for the free electromagnetic field, whereas the zeroth order 
term is given by 

TZM\ = -d»A d»A»+ l -d,A u d»A» 



i (V^°) 2 = ^aWM (32) 



2 

and represents the energy density of the periodic background field. Due to the 
projection to a particular sector of the theory, the additional term 

T A 00 [A,c] = (c - A»(x)A,(x))\ 

'(IT 1 I duA^u)A^u) - A^{x)A fl {x)]x 

= -a 2 (l-2cos 2 (£x 1 ))X (33) 

occurs. It is easy to see that this term does not contribute to the energy, but 
gives a non-vanishing periodic contribution to the energy density. Finally, the 
term 

T oo = _ l _ d p(^ J o^ +1 p rr og D ^ A + a ^ (34) 

represents the contribution of the fermions to the energy density. The first 
term on the r.h.s. gives vanishing contribution to the total energy if charge 
conservation is required, whereas the second term accounts for the kinetic 
energy of the fermion system. At one-loop order, we have to substitute a M = 
in the 'Hamilton operator', i.e. write 

H D (A) = -i 7 V^-7V^W+7 m. (35) 

Since Aj(x) = in our case, we obtain Hu{A) = —vf^&j + 7°m, i.e. the 
free 'Hamilton operator'. 

The vacuum expectation value of T 00 in the presence of the background 
field A^ is defined as 



(T 00 (x)) = fd\f VaVipViljT 00 {x)e i{Sx+SD) 1 1 dX f VaV^V^e i{s * +SD) 

(36) 

Similarly, the vacuum expectation value of Tjy ee in the absence of the back- 
ground field (Tj?° ee (x))o is given by substituting A M (x) = in Eq.(^) and 
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replacing T 00 (x) by T^ ee (x). Then, the Casimir energy of the vacuum due to 
the background field A^(x) ^ is 

E c = Jdx 1 (T 00 (x 1 ))-Jdx 1 (T^. ee (x 1 ))o 

= J dxi [(Z V/2) sm 2 {l Xl ) + (f T °#oo*> - (% iJ DO *)o] (37) 

This expression of the Casimir-energy reminds one on the expression for the 
energy of a system of electric charges in classical electrodynamics. There, the 
energy is the sum of the energy of the electromagnetic field and the kinetic 
energy of the charges [|J. In our case the Casimir energy is the sum of the 
background electric field and the change of the relativistic kinetic energy (in- 
cluding the rest mass) of the Dirac-sea due to the presence of the background 
field. 

Eq.([3^) can be rewritten by the help of the generating functional 

Z'qed[Ac,J,V,v] as 

00/ \\ _ ft Iryl r T n l\ rpOO 



(T 00 (*)> = {l/Z f QED [A,c })T™(x)Z' QED [A,c,j,r),fj} (38) 

where T^(x) denotes the operator obtained from T 00 by replacing the fields 
a^(x), ip{x), and ip(x) by the operators 8/6j^(x), 6/8fj(x), and —S/Srj(x) and 
the variable A by iQ~ 1 d/dc. 
Then, we find 

(T A 00 > = -a 2 e 2 (2cos 2 (ex 1 )-l). (39) 

Furthermore, the expectation value {T^ 2 [°\) * s equal to the energy density of 
the free electromagnetic field. Indeed, it holds for the second derivatives 

8 2 Z' em /SUy)SUx)\ = iDr(x,y)Z' em \ J=0co (40) 

where D^ v {x,y) tends to the free propagator D^ u {x,y) in the infinite volume 
limit, since the last term on the r.h.s. of Eq. (|T8| ) is of the order Qr x . Conse- 
quently, the pure electromagnetic contribution to the Casimir energy density 
is given by 



(*) = «,o^]+^ 00 [Ac ]) 

1 



^00 r T\ _i T^OOr 
em,ol/\l "r 

,2«2 „„„2/ 



a 2 ^cos 2 (^i). (41) 
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It is more cumbersome to evaluate the fermionic contribution cf{x) = 
(Tp°) — (Tjy ee )o to the Casimir energy density, where the energy density of the 
free Dirac vacuum is subtracted. Since QED in dimensions 1 + 1 is superrenor- 
malizable, the difference turns out to be UV finite without further renormaliza- 
tions. We perform the evaluation of the fermionic part of the Casimir energy 
in the second quantized formalism. Then we can write ep = (0| : T F ° : |0) 
where |0) is the 'interacting' vacuum (in the presence of the periodic back- 
ground field) and : . . . : denotes normal ordering with respect to the normal 
vacuum (in the absence of the background field). The evaluation is performed 
in the following steps: 

1. The fermion field ip is expanded in terms of the eigenspinors f^ ks \xi) 
and g( ks \xi) of the Dirac Hamiltonian Hd(A) belonging to the energy 
eigenvalues eks and — e^s, resp. Here the quasi-momentum k G \—&/2, £/2] 
and the integer s > enumerating the bands are introduced. 

2. Then the creation and annihilation operators a ks , b ks and a& s , bk s of these 
stationary single particle states are expressed as linear combinations of 
the creation and annihilation operators A ks) B ks and A ks , B ks of the free 
fermion states F^ kr \x\) and G^ kr \x\) of energies e k J = Jm 2 + (k + £r) 2 
and — e kr , resp. In terms of the latter, the normal ordering is performed. 

3. Finally, the normal ordered operator is reexpressed in terms of the oper- 
ators a ks: b ks and aks, b ks and the vacuum expectation value with respect 
to the vacuum in the presence of the background field is taken. 

Thus, one arrives to the following expression 

Tf °(x) = e F (x)+^{al,(t)a p (t)B p ;-Bi + bl(t)bAt)A p ;-A p _ 

pp' 

+b p ,(t)a p (t)A p ; ■ B p _ + a^{t)b\{t)Bi ■ A p _) (42) 

where, with p = (ps), p' = (p's'), 



and 



Ai = &Y + al p \ 


K=py+Pi p \ 


A p _ = a F — a p Gl 


Bl=l%-I%, 


^2e kr a p _ ksr F kr {x), 


ps \~^ pksr /o 

a G =2s 6 kr a + G 



(43) 



ps V~~* Pksr t-i / \ ps \~~\ pksr /^ i \ 

«f =L £ fcr«- F kr {x), a' G =2_^e kr a F + G kr (x), 

kr kr 
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0? = 


EV- fa %(4 




kr 


~* ps 


= Y,a*J ksr FZ r (x), 



kr 



a 



kr 
kr 

with the constant coefficients 

a p _ ksr = fdxF* kr (x)g ps (x), 
(3 pksr = fdxF* kr (x)f ps (x), 



kr 



i\ 



ps 



kr 



(44) 



a p + ksr = dxG* kr (x)g ps (x), 
Pl ksr = fdxG* kr (x)f ps (x). 



(45) 



The latter are the overlap integrals of the eigenspinors in the presence of the 
background and those in the absence of the background. Furthermore, the 
time dependent creation-annihilation operators are introduced: 



, ps\J'j — O'ps^ > 

The constant term 



™ps\£) — "p S € 



lErjst 



~* ps ps -■ 

a G ■ a F — a 



ps ps 



a 



G 



(46) 



(47) 



ps 



represents the Casimir energy density of the fermion vacuum. The last two 
terms of Eq. (|4"7| ) cancel. 

Performing the calculation results in: 



e F (x) 



E e\ r ^r r) +n krr) )oos{{r-r'yx 



ke[-±,±]rr>eN 

,0 \ r-ikrr') 



where we introduced the following notations: 

j-(krr') = y^ jj{-kr')jj(-kr)jj{-kr') v {-r'ks)jj(-kr) v (-rks) 
sez 

j-(krr') = v^ y(- kr )y(-kr')y(-kr) u (-rks)y(- kr >) u (-r'ks) 

s&Z 
j-.(krr>) _ y> y(-fcr , ) V -(fcr)yr(-fcr') w (-r'fc») v -(fcr) t; (rfc«)_ 



(4J 



(49) 



s£Z 



[/( fcr ) and V^ denotes the eigenspinors of the free Dirac-equation for the 
positive and negative energy eigenvalues, respectively. It is straightforward to 
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Figure 4: Casimir energy per unit volume of the periodic vacuum in QED 2 



establish from Eq. 
negative: 

E, 



that the volume integral of the energy density is not 



krs 



-kr 



\u 



(—kr) (rks)\2 



P U P 



+ W f 



(—kr) (rks)\2 



P 



'V 



ft 



> 0. 



(50) 



By numerical calculations we were convinced (see Fig. |j) that this expres- 
sion of the energy only vanishes in case of a=£=0. This means that the Casimir 
energy E c of the vacuum is always positive if a non-vanishing periodic electric 
background field is assumed. Thus, the vacuum of QED 2 does not favour a 
periodic mean field energetically with respect to the normal vacuum. 



6 Necessary condition of energy minimum 

The vacuum expectation value A°(x) is defined by the minimum of the energy 
functional 



TE[A°(x)} = Jdx(T 00 (x)) = Jdx((TZ(x)+T™(x)+T J 



° F °(x 



(51) 
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The second term vanishes due to the explicit value of Co- Thus, the necessary 
condition of the energy minimum takes the form 



6 



5A°(x) 



dy(TZ,olA°(y)] + (T F °(y)))=o. 



The functional derivative of the first term gives 
S 



5A°(x) 



dy(TZA AU (y)] 



V z A u (x) 



(52) 



(53) 



At one-loop order the fermionic contribution to the energy can be rewritten as 

te f = /dy(r°°( y ); 

= Z?[A, 0, 0] J dyT™ p (y) Z F [A, 77, fj] 

Taking its functional derivative we find 

5 In Z F [A, 0,0], 



r)=ri=0 



(54) 



5__ 

5A°{x 



-TE, 



SA°(x) 



'-TE, 



+Z^L4,0,0] UyTfiy) 



5Z F [A,r],T]} 



5A°(x) 



rj=r)=0 



i/dyfkfT^T 



s 



2 y \Sr]{y) 5i](y) 



8rj(y) 



DO 



1 U H, 



5 



8rj(y) 



n 



5G 



5A°(x 



-77 e 



i(r)Gr)) 



-itr (y°G(x,x)) + c°(x) 



(55) 



with 



c (x) = - dy[tT(7 G(y,xh D° y G(x,y)) 



tr( 7 G(x,y) 7 WG(?/,:r))]. 



(56) 



Thus, we find the following equation for the field A°(x) minimizing the energy: 

d v F 0u = U(x)~f°ij(x)) + c (x). (57) 



The first term on r.h.s. is the expectation value of the charge density j° = 
—itr('j G(x,x)), determined via the propagator G(x,x) as a given functional 
of A°(x), and a similar statement holds for c°(x). 
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A similar equation appears in the classical case when a certain charge 
distribution moves in an electromagnetic field. Integrating out the effect of 
the charged particle distribution it will result in a polarisation charge density 
term beside the common charge density in the equation of motion for the 
scalar potential. In Eq. ( |57D the term c°(x) corresponds to a polarisation 
charge density, so minimizing the energy with respect to A°(x) leads to an 
equation similar to that of Poisson's equation in a polarised medium. Thus 
Eq. ( |5?D can be solved directly for A (x) in principle without the need for 
solving any other equations. On the contrary to this the mean field equation 
fl2T|) obtained by extremizing the effective action does not take the polarisation 
of the vacuum into account, in order to do this we have to solve a system of 
operator equations as well. 

7 Conclusions 

The energy density of QED2 in the presence of a periodic mean field is de- 
termined from the energy- momentum tensor. For this purpose a projection 
method is worked out which is applied to treat the mean electromagnetic field 
self-consistently. The projected effective action and the energy density of the 
vacuum are derived at one-loop order, whereas the interaction of the electron- 
positron field with the periodic mean field is treated exactly. It is established, 
that the negative of the effective action must not be regarded as the energy 
of the system if the background field is not constant. It was shown that the 
necessary conditions of the extremum of the effective action and the minimum 
of the energy lead to different equations for the vacuum expectation value of 
A . Eq. ( p7|) obtained by extremizing the effective action is the vacuum expec- 
tation value of the Poisson's equation that does not include the polarisation 
of the vacuum due to one-loop radiation corrections. Those are accounted for 
by separate operator equations for the quantum fluctuations of the electro- 
magnetic field and for the Dirac field. On the other hand Eq. (|57|) obtained 
by minimizing the energy functional includes the polarisation effects of the 
vacuum. 

The expectation value of the component T 00 of the energy-momentum ten- 
sor was determined as the function of the amplitude a and wave number / of 
the static periodic scalar potential A°(xi) = acos(lxi). It is found that the 
vacuum configuration with this periodic electric mean field is not favoured en- 
ergetically compared to the normal vacuum. The volume integral of the energy 
density plotted against the amplitude a in Fig. f| shows that the energy of the 
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system increases with increasing a monotonically. 

The result obtained contradicts to our naive expectation discussed in the 
Introduction. Possibly, the reason is that the naive picture mentioned there 
does not take into account the uncertainty principle i. e. for a certain I the 
energy spread of a wave packet localized in an interval of ~ l/£ could be much 
larger than the amplitude of the potential and of course in this case our naive 
picture is not valid any more. We set now forth our work looking for periodic 
ground states at finite chemical potential. 

A Solution of the Dirac equation in periodic 
external field 

A.l. Relativist ic Bloch waves 

To get the projected effective action at one-loop order in Eq. (^) we must find 
the fermionic single-particle energies, so we have to solve the Dirac equation 
in the presence of the potential flTT|): 

(i-f{d fl -ieA fl )-m)^ = 0. (58) 

We look for the solution of Eq. ( |58|) in the form of Bloch waves corresponding 
to the energy eigenvalues E: 

,/. _ p -ie P st ipx V^ inlx _ I e PS J { x )i ^ E = e ps > U , , 

Wa e e y Uane - \ e **.y»( x ), if E = -e ps > { } 



Inserting fl59D into (|5q) we find 



£ ( Man (e 7 -(p-n/)7 1 -ml) + ^° (u an+1 +u an ^)y^ +nl >- iEt = 0. (60) 

n=— oo ^ ' 

We get non-trivial solutions when the determinant of the matrix appearing next 
to the Dirac-spinors equals zero. If we had solved the Schrodinger equation 
in the presence of this sinusoidal potential, the form of the matrix would have 
been tridiagonal which means that the diagonal and the neighouring diagonal 
elements are nonzero. In relativistic case the matrix elements are replaced by 
d dimensional 7 matrices but the structure of the matrix remains unchanged. 
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A. 2. Eigenvalues, eigenspinors 

We cannot get the eigenvalues analitically because of the complicated structure 
of the matrix but we can determine them as precisely as we wish. In numerical 
calculations we work with matrices with finite dimension. Using the well- 
known identity || 

detL4] = det (^ ®) = det[P] det[5 - RP^Q] (61) 

we can reduce the problem of calculating the determinants of these d(2n + 
1) x d(2n + 1) dimensional matrices to calculating four dimensional matrices 
by identifying the upper left matrix element of the matrix under investigation 
with P in Eq. (|6T|). (n denotes the number of u a ^-s taken into account in Eq. 
(|60D.) Using identity (|6l| ) In times we get a product of 2n determinants of 
four dimensional matrices. We computed the determinant of Eq.(^) as the 
function of E and determined its zeros, corresponding to the energy eigenvalues 
E = ±e 

In order to evaluate — r we have to sum the eigenvalues — e ps < and 
extract from it the sum of the negative eigenvalues of the free Dirac equation. 
This difference will depend on the accuracy of the eigenvalues, the number of 
the eigenvalues taken into account in the sum, the size of the chosen matrix 
and, of course, the parameters of the potential, the amplitude a and the wave 
number I. We have to be convinced of the stability of the numerical calculation. 
We increased the number of members in the sum and the numerical accuracy 
of the determination of eigenspinors by choosing larger matrices as far as we 
have seen that the energy difference does not change significantly. 

To get the expectation value of the component T 00 of the energy-momentum 
tensor we also have to determine the eigenspinors of the Dirac equation (p8|) . 
We calculated them with the help of the eigenvalues by solving a system of 
homogeneous linear equations (|60| ) for the eigenspinors. 
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